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We explore systems of pulse-coupled oscillators beyond the mean-field limit [R.E. Mirollo and S.H. 
Strogatz, S1AM J. Appl. Math. 50, f645 (f990)] by means of a manageable description which 
leads to a great simplification of the dynamics. Sufficient conditions for synchronization are exactly 
obtained for a ring with nearest-neighbor directed interactions, turning out to be the same as in the 
all-to-all case. Our analysis shows new synchronization mechanisms occurring when the interactions 
are local. 



PACS numbers: 05.45.+b, 87.10. +c 

Nonequilibrium spatially extended dynamical systems, 
consisting of many degrees of freedom, constitute one of 
the most frequent manifestations in nature, which how- 
ever still require new approaches to be understood [Q. 
Very remarkable in them is the spontaneous emergence 
of order, i.e., self-organization, by means of local, short- 
range interactions that cooperate to build correlations 
over all spatial and temporal scales || ■ Among all forms 
of order in space and time, synchronization and periodic- 
ity are the strongest ones, respectively, and play a crucial 
role in many physical, chemical, and biological systems 



Systems behaving in this way are in fact populations 
of oscillators coupled between them; in particular, out- 
standing examples make it clear that in most living sys- 
tems the interaction between units cannot be viewed 
as continuous in time but rather as discontinuous and 
episodic, in the form of sharp pulses from one oscillator 
to another. Pacemaker cells in the heart [||, neurons in 
the visual cortex Q , and fireflies that flash on and off si- 
multaneously in periodic intermittency || provide strik- 
ing illustrations of spontaneous synchronous evolution in 
nature, all of them with pulsatile interactions. 

The discontinuous character of the coupling makes 
remarkably difficult the mathematical analysis of these 
systems. In the case of a network of all-to-all pulse- 
coupled oscillators, Mirollo and Strogatz, however, were 
able to overcome the hardness of the problem establish- 
ing a sufficient condition for synchronization ||. Their 
result is based on the idea of absorption, which means 
that once two oscillators get synchronized they never 
desynchronize. In this way, they proved that, no mat- 
ter the initial conditions, synchronization arises cooper- 
atively by a sequence of absorptions, in a self-organized 
process where no cell leads the population but it is the 
whole that coordinates the overall activity to achieve self- 
synchronization. 

In contrast to the Mirollo and Strogatz 's scenario and 
other situations with all-to-all, or mean-field, interac- 
tions Jl0|-|l3t, for more realistic networks jL4j, in which 
the firing of one unit does not modify the state of all the 



others in the same way, the concept of absorption is no 
longer valid and the synchrony achieved by a subset of 
oscillators can be lost p^pq j. Thereby, in these cases 
one would expect a weaker tendency to synchronize, as 
occurs in other large networks of oscillators. Numeri- 
cal studies O], however, indicate that the conditions for 
synchronization in networks with reduced connectivity 
seem to be the same as in the all-to-all case. Despite 
this counterintuitive result was already conjectured by 
Mirollo and Strogatz ||, up to now there is no proof of 
synchronization in any model of pulse-coupled oscillators 
beyond the mean-field assumption. 

In this Letter we address precisely the problem of es- 
tablishing sufficient conditions for synchronization to oc- 
cur in networks of pulse-coupled oscillators with local 
connectivity. In this regard, we have exactly derived 
them for the Mirollo and Strogatz 's model with nearest- 
neighbor directed coupling in one dimension, turning out 
to be the same as in the mean-field case |)| . This kind of 
feed-forward connectivities arise frequently in the model- 
ing of some layered structures of the brain |Q , whereas 
in loop- like circuits, they are relevant for cardiac arrhyth- 
mia and the brain clock Jl9| . 

Consider a population of N relaxation or integrate- 
and-fire oscillators, where the units evolve in time inte- 
grating an excitatory driving signal until the state vari- 
able of some of them reaches a threshold, firing an instan- 
taneous zero-width pulse to their neighbors, and then 
relaxing to a lower-value state, where the driving acts 
again. In general, any system of oscillators can be de- 
scribed in terms of the phases fa of its units |]l3| , f20f . For 
relaxation oscillators fa measures the normalized elapsed 
time since the last firing of unit i if there had been no in- 
teractions. The time evolution of the units is then given 
by 

if^<IVz^^ = i, i = l,2,.. .N, 



(1) 



if fa > 1 



4>n(i) -> fai(i) + A(0„(i)) if > 0, 
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where T is the period of the oscillators and A(<j>) the so- 
called phase response curve, which gives the phase shift 
of a unit when it gets a pulse from a neighbor. Neighbors 
of unit i are denoted by n(i) and for a one-dimensional 
directed lattice with coupling "from left to right" , n(i) = 
i+1. Periodic boundary conditions are taken assuming 
N + 1 = 1 and ee N. 

Notice that T and A(</>) are the same for all the units, 
which therefore arc identical and identically coupled. Ob- 
serve also that a unit at the reset point, <f> = 0, does not 
feel the incoming pulse; this is a way to include the ef- 
fects of absolute refractoriness. In addition, as in the 
Mirollo-Strogatz model, we only deal with excitatory in- 
teractions, i.e., A(0) > 0, V0; this makes the interaction 
process evolve in the form of avalanches. 

We will show that a sufficient condition for com- 
plete synchronization is obtained when the phase re- 
sponse curve is an increasing function of the phase, i.e., 
A'((^>) > O,V0. Our proof of synchrony consists of two 
differentiated parts. The first one establishes that once 
an oscillator has fired triggered by the firing of its neigh- 
bor, both units will fire at the same time forever, even 
though they may lose their mutual synchrony between 
firings. We call such a process a capture. The second 
part shows how for almost all initial conditions there is a 
capture for each pair of neighboring units. In this way, as 
eventually every pair of neighbors are mutually captured, 
the firing of one unit will ensure the firing of the rest and 
the complete synchrony of the whole population will be 
achieved, with all the units evolving with the same phase. 

In order to prove our first claim, let us consider two 
neighboring units, labeled by i and i + 1, which have fired 
at the same time. Since both units have the same free 
time evolution, they can only desynchronize when the ith 
one gets a pulse; when this happens, at time a, its phase 
<\>1~ is instantaneously increased by an amount A(</>" - ). 
Here the superscript of the phase indicates time. Notice 
that just before the pulse both phases are equal, i.e., 

<t>r = c+i • (2) 

In contrast, just after the pulse the phase of the (i + l)th 
unit remains at the same value, in principle, whereas the 
phase of the ith unit changes to 

#+ = #-+A(#")- (3) 

At this point, two different situations may occur depend- 
ing on whether the ith unit crosses the threshold as a 
consequence of the received pulse or not. 

If <j>1~ + A(0?~) > 1 the unit crosses the threshold and 
fires, increasing then the phase of its right neighbor by an 
amount A(^ 1 ). Since A(^) = A(0?~), the (i + l)th 
unit also crosses the threshold, firing at the same time as 
the ith one. 

If 4>1~ + A(0? _ ) < 1 the ith unit docs not reach the 
threshold. Therefore, after receiving the pulse its phase 



increases continuously in time until either it reaches the 
threshold or it gets another pulse. In the former case, the 
unit fires at time b = a + T[l- - A(0? - )], sending 
a pulse of strength A(^^ 1 ) to i + 1. Notice that as 

^- 1 = i-A(0n, (4) 

the phase of i + 1 after the pulse will be 

^^l-A^n + A^x). (5) 

Assuming that A(cf>) is an increasing function of its argu- 
ment, 4>i+i is greater than 1. Consequently the (i + l)th 
unit crosses the threshold, firing then at the same time 
as its left neighbor. 

Let us now discuss the case in which the ith unit re- 
ceives two pulses. Indeed, as the strength of the pulses 
is not fixed, a unit can overtake its right neighbor; this 
is what makes a unit get more than one pulse while its 
phase is evolving between and 1. For the sake of gen- 
erality we also consider that a set of connected units on 
the left of i also get two pulses, except the leftmost unit 
of this set, denoted by L, which only gets a single pulse. 
The existence of a unit receiving only a single pulse is 
guaranteed by the periodic conditions and by the fact 
that at least i + 2 gets only one pulse before i + 1 can 
reach the threshold. After sending the first pulse to unit 
j, j — 1 is delayed with respect to j and consequently 

< f/- 1 , Vj s.t. L < j < i , (6) 

where the superscript (Zj-i indicates any time between 
the first firing of the (j — l)th unit and the time in which 
j — 1 gets its first pulse in this cycle. In order for a 
unit j to get the second pulse, the second firing of its 
left neighbor must occur before the second firing of j. 
The sequence of firings of this set is then from left to 
right, starting on the Lth unit. Therefore, just before 
and after the Lth unit gets its only pulse at time a, its 
phase must be smaller and greater, respectively, than the 
corresponding ones of the remaining units in the group, 
i.e., 

4>r < 4>T < 4>l~ + A(0D , Vj s.t. l < j < i + 1 . 

(7) 

Once the Lth unit reaches the threshold at time b = 
a + T[l - ^ - A(0° - )], all units from L to i + 1 fire 
because 

$ + = 1 + 4>T 4>T + A[4> b f] A[#"] > i . 

Vj s.t. L < j < i + 1 , (8) 

since A'O) > 0, V0 and <j> b r > <f~ , Vj s.t. L < j < 
i + 1. 

Our analysis then indicates that once two neighboring 
units fire at the same time, they always remain firing 
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FIG. 1. Spatiotemporal evolution of a ring of 100 oscilla- 
tors for A(0) = e + with e = 0.01 and 7 = 0.3. Black and 
white colors stand for minimum and maximum values of the 
phase <pi, respectively. Notice that the coupling is "from left 
to right" . 

at the same time if the left unit only gets one or two 
pulses. The case in which the left unit gets three or more 
consecutive pulses is not allowed in this process because 
a unit reaches the threshold and fires when it gets the 
second pulse, as we have shown. 

Concerning our second statement, we now demonstrate 
that the set A of initial conditions of an arbitrary oscil- 
lator for which it will never fire at the same time as its 
left neighbor has Lebesgue measure zero. To be precise, 
for a unit i + 1 and given an initial configuration of the 
rest of the population, the set A® +1 is defined as 

A° i+1 = {$+1 s.t. i will not capture i + 1 and 
i will not be captured in an avalanche triggered by i + 1, 
provided that the rest of the population 
is in an initial state {(^j-^i+i }. 

In particular, if = 1, A° i+1 C s.t. < < 

1 - AO*)}, where </>* is defined by + A(<j)*) = 1; 
otherwise, i captures i + 1 immediately. In essence, the 
proof shows that for almost all the initial conditions of a 
unit, and irrespective of the state of the rest, it cannot be 
confined forever in the interval (0, 1 — A(^>*)) receiving 
an infinite number of pulses from its neighbor. 




FIG. 2. (a) Same situation as in Fig. hi but e = 0.1 and 
7 = 0. (b) Same situation as in Fig. hTbut e = 0.7 and 
7 = -0.1. 

Due to the characteristics of the dynamics, the time 
evolution of any unit can be decomposed into two parts. 
On the one hand, the free time evolution in which the 
unit does not get any pulse is given by 

4>(t') = **'-*[#*)] = U(t) + ^) mod 1 , (9) 

where t and t' are times between two consecutive pulses. 
Notice that time is now indicated by the argument of the 
phase. On the other hand, the dynamics of the unit is 
coupled to its neighbors. Thereby, when its left neighbor 
fires the unit gets a pulse and increases its phase by an 
amount A(0), provided that it is neither brought to the 
threshold nor refractory, which, by definition, is always 
fulfilled for the time evolution of the set A . Then, the 
new phase just after the firing is 

r +1 (t) = hint)} = r(t) + A[«r(t)] , (10) 

with the superscript n indicating now the number of 
pulses received by the unit. 

It can be seen that the measure [i of the time evolu- 
tion of the set A remains invariant during the free time 
evolution by using that 

dFt '-y =l for o< 0<1 . (id 

d<p 

Indeed, if A{t') denotes the state of the set .4° after a 
time t' , 

f,[A(t')\ = [ dpW)] = f dFt ~Jf {t)] d^(t)} 

JA(t') JA{t) d<f>{t) 

= [ d^{t)\ = ■ 

JA(t) 

In contrast, when the unit gets a pulse, the measure 
evolves in such a way that the inequality 
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(12) 



holds, with A' = min[A' (</))] and A n+1 being the result- 
ing of the evolution of the set .4° after the n + 1 pulses. 
This follows directly from both 



n+l\ 



dh((j) n ) 



dnW) (13) 



A" 



and dh(cj) n )/d(j) n = 1 + A'(<p n ) > 1 + A'. If A(</>) is 
an increasing function of its argument, its derivative is 
positive and then 1+A' is greater than 1. In consequence, 
if the unit always gets the pulse at <f> G (0, 1 — A(cj)*)), 
when n goes to infinity the inequality 



(i(A 



(1 + A')" +1 



(14) 



implies that the set .4° has zero measure since /i(A n+ ) 
is bounded by 1 — A((f>*). 

In other words, the flow in phase space remains incom- 
pressible during the driving; however, if A'((j>) > 0, V0, 
the coupling acts in a more complicated way, expanding 
the volume elements with each firing, except in the case 
of a capture, where the volume contracts. This is the pro- 
cess which eventually dominates, leading the population 
to the complete-synchronization attractor. In addition, 
it is easy to show that this condition is fully equivalent 
to the one obtained by Mirollo and Strogatz |2(| . 

To illustrate the essentials of the emergence of syn- 
chrony in systems with explicit spatial dependence, in 
Fig. |l| we have depicted the results obtained from nu- 
merical simulations for a directed ring when the condition 
A'(<p) > 0, V</>, is satisfied. This figure shows the time 
evolution towards the synchronized attractor, starting 
from a random initial condition. It is evidenced how the 
complete synchrony of the whole population is achieved 
through a coarsening process in which the domains of 
captured oscillators grow until eventually all the units 
are evolving with the same phase. In contrast, in Fig. 
U we show the time evolution for two situations wherein 
synchronization does not appear. Fig. ||a corresponds to 
a case in which A'((f>) = 0, V</>. In this situation, starting 
from a random initial condition, perpetual disorganiza- 
tion occurs; the system remains evolving periodically in 
time in a spatially disordered state. The situation shown 
in Fig. ||b corresponds to a typical pattern reached after 
a long transient (not displayed) for A'(cj)) < 0, V0. In 
this case, synchronization does not occur although the 
final pattern looks more regular than that of Fig. [2^,. 

In summary, we have established, for the first time, suf- 
ficient conditions for synchronization in a general model 
of pulse-coupled oscillators with local interactions, en- 
suring that, on a directed ring, they are the same as 
in the mean field case. The simple model of rhythmic 
oscillations in biological systems we have presented il- 
lustrates the fundamentals of emergence of organization 



in nonequilibrium complex systems led by local rules, in 
contrast to previous work devoted only to global cou- 
pling. Our analysis, then, elucidates new mechanisms 
taking place to attain synchrony when the spatial struc- 
ture is explicitly considered 
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